Four-level and two-qubit systems, sub-algebras, and unitary integration 
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Four-level systems in quantum optics, and for representing two qubits in quantum comput- 
ing, are difficult to solve for general time-dependent Hamiltonians. A systematic procedure is 
presented which combines analytical handling of the algebraic operator aspects with simple solu- 
tions of classical, first-order differential equations. In particular, by exploiting su{2) © su{2) and 
sii(2) ffi su(2) ® u{l) sub-algebras of the full 517(4) dynamical group of the system, the non-trivial 
part of the final calculation is reduced to a single Riccati (first order, quadratically nonlinear) equa- 
tion, itself simply solved. Examples are provided of two-qubit problems from the recent literature, 
including implementation of two-qubit gates with Josephson junctions. 

PACS numbers: 03.65.Yz, 05.30.-d, 42.50.Lc 



Four-level systems are often of interest in quantum op- 
tics and in multiphoton processes. As examples of re- 
cent literature, see, for instance. Added interest is 
provided by recent applications for quantum information 
and computing where two qubits constitute a four-level 
system. Indeed, two qubits and associated logic gates 
have been practically realized through Josephson junc- 
tions 3], the coupling Haniiltonian viewed as a 4 x 4 
matrix. 

In recent papers, we have developed a "unitary inte- 
gration" technique for solving time-dependent operator 
equations such as the Schrodinger or Liouville equation 
for the state or density matrix of a quantum system . 
The method also extends to more general master equa- 
tions when dissipation and decoherence are present j^. 
The evolution operator is written in the form of a prod- 
uct of exponentials, each exponent itself a product of 
a classical function of time and one quantum operator. 
When these operators form a closed algebra under com- 
mutations, first-order time-dependent equations can be 
derived for the introduced classical functions, and their 
resulting solutions allow the construction of a complete 
solution of the quantal problem. The method works best, 
of course, when the number of operators in the algebra 
is small. For one qubit, where the most general Hamilto- 
nian involves only four operators, the three Pauli matri- 
ces and the unit matrix, the su{2) algebra has just three 
elements. Upon including dissipation and decoherence, a 
3x3 matrix for the elements of the density matrix in- 
troduces, at its most general, eight elements of an su(3) 
algebra |5j- For special situations when an su{2) sub- 
algebra suffices, such decoherence problems can also be 
treated through just three classical functions whose defin- 
ing equations further reduce to just one non-trivial Ric- 
cati (first-order, quadratically nonlinear) equation (or, 
equivalently, to a linear second-order differential equa- 
tion). 

The aim of this paper is to examine two-qubit or four- 
level problems analogously. The complete description in 
terms of the 15 operators of the symmetry group SU (4) 



(together with the unit operator) requires 15 exponen- 
tial factors in the evolution operator which is cumber- 
some. However, judicious choices of sub-algebras reduce 
the problems considerably. These may be of interest since 
they describe several of the two-qubit logic gates and 
four-level systems in the recent literature. Specifically, 
we present here the cases pertaining to two sub-algebras, 
su{2)(Bsu{2) and su(2)©su(2)©m(1). In both instances, 
the same single Riccati equation for a single function is 
all that remains to be solved (through, for example, a 
popular Mathematica package 0) to provide complete 
descriptions of such systems. Other sub-algebras, namely 
so(5) with ten operators and su(3) with eight, will be the 
subject of future work. 



ALGEBRA OF OPERATORS FOR A 
FOUR-STATE SYSTEM 

SU{2) X SU{2) sub-groups 

The sixteen linearly independent operators of a four- 
state system can be chosen in a variety of matrix repre- 
sentations. One choice, familiar in nuclear and particle 
physics applications are the 15 matrices of su(4). 
This set is suitable for examining su(3) and su(2) sub- 
algebras involving 8 and 3 operators, respectively, but 
for our interest in this paper, an alternative representa- 
tion proves more convenient. This starts from two qubits 
and employs their individual Pauli matrices, 6 in number, 
along with tensor products, 9 in number: ^a, , ^ct® ^t. 
Together with the unit 4x4 matrix, these 16 matrices 
called Oi are explicitly tabulated in an earlier paper 
(A slightly different arrangement of the same operators, 
called Xi , with much useful group-theoretic discussion, is 
in 9] ) . They have the advantage for the unitary integra- 
tion procedure (to be described below) that the square 
of each is, to within a factor, the unit matrix. Expo- 
nentiation of these operators is thereby rendered simple, 
involving just two terms. 
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TABLE I: Table of commutators. With operators Oi in the first column and Oj in the top row, each entry provides the 
commutator [Oi ,Oj]. 



Table I gives the commutators [Oj, Oj] and is useful for 
picking out sub-algebras. For example, [02,05] = iOg. 
The four states in this representation may be described 
either in the four-state language as the column vectors 
(1,0, 0, 0), (0, 1, 0, 0), (0, 0, 1, 0), (0, 0, 0, 1), or in spin no- 
tation as I TT), I Ti), I it), I ii), respectively. 

There is a large set of six-dimensional sub-algebras 
formed out of two triplets obeying su{2) or, equiva- 
lently, so(3) angular momentum commutation relations. 
One such set, (O2, O5, Oe; O3, O9, Oio), is the trivial 
su{2)(Bsu{2) subalgebra of local operations on two qubits 
a and f . There are many equivalent sets generated by a 
Similarity Transformation (ST) acting on the elements of 
a subalgebra by right-left multiplication by an arbitrary 
element of the full S'C/(4) group. Let be a unitary 
4x4 matrix, and x an sm(4) algebra element. Then, ST 
acts as 

x' = WxW^ = WxWK 

In this section we will use only a discrete subgroup of 
all possible ST transformations, which have the property 
that the set of 15 operators Oi is transformed into the 
same set Oj, i.e. the action of ST is equivalent to a 
permutation among the Oi. In this case the su{2)®su{2) 
subalgebras, obtained by ST transformation of a, r can 
be identified by direct analysis of subgroup structures 
using Table 1. 

Subgroups different from the trivial a, r one can 
be generated only by two-qubit entangling operations, 
such as transformation to the Bell basis (or magic 
Bell basis ^3)- One can easily verify that the set 
(O9, Os, O14; Oe, On, O13) is such a subalgebra. That is, 
the operator sets (xr, GyT^, OyTy) and (ctj,, UzTx, OxT^) are 
two mutually commuting triplets with su{2) commuta- 
tion relations. The ST transformation with the unitary 
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passes from the a set to {Tx,ayTz,ayTy) and from t to 
{ay,azTx,(JxTx)- The basis corresponding to Eq. Q is 
the Bell basis, (1,0,0,1) = | TT> + I ii>, (0, 1, 1, 0) ^ 

I Ti ) + | it), (0,1,-1,0) = I ti)- I it), (1,0,0,-1) = 

I TT) ~ I ii). which are simultaneous eigenstates of axTx 
and (TyTy [lO|. 

All such subalgebras are maximal su{2) © su{2) sub- 
algebras of su(4), in the sense that, if any operator from 
Table I is added to the six already chosen, the algebra 
closes (as a Lie algebra) only in the full sm(4) algebra. 
Another perspective on the above Bell basis is provided 
by considering a four-level system with nearest neighbor 
couplings as in the Hamiltonian, 



H = 



/ a \ 

a P 

/3 7 

V 7 / 



(2) 



Such a matrix may be viewed as built up of three sym- 
metric 4x4 matrices rrijk, with unit entries in row (col- 
umn) j and column(row) k and zero everywhere else, so 
that H — ami2 + /3to23 + 7'Ti34- Starting with these 
three real, symmetric matrices, we define two non-zero 
commutators between them, TO13 = —«[mi2, 77123], W24 = 
— i[?Ti23, W34] to give two imaginary, anti-symmetric ma- 
trices, and a final real symmetric 77714 = ^7[77ii3, 77134] 
to complete a full 6-dimensional so(4) algebra. Using 
the well-known relation between 5*0(4) and SU{2) x 
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SU{2) groups, this sub-algebra decomposes into two 
sets of mutually commuting triplets, mi2 + ?ti34,toi3 — 
"i24, rni2 - and TO12 - 77134, 7*7,13 + "724, "^23 - "^14- 
These linear combinations correspond to the subalgebra 
(O9, Os, O14; Oe, On, O13) of the full sm(4) algebra. Such 
a construction generalizes to 77-level systems, 77 > 4, with 
nearest neighbor couplings. The matrices 777^^ defined as 



{rrijk] 



pq 



\j-k\- 



-1) 



k-j- 



provide 77(77 — l)/2 generators of a 50(77) algebra. These 
form a sub- algebra of the 77^ — 1-dimensional su{n) alge- 
bra although, for 77 > 4, the 50(77) no longer reduces to a 
product of S7i(2) algebras. 



5'I7(2) X 5'I7(2) X ;7(1) sub-groups 

A class of maximal seven-dimensional sub-algebras 
of 57t(4), different from the one considered in the pre- 
vious subsection, can be constructed by taking linear 
combinations of 0„ generators. Inspection of Table 
I shows that the first seven elements form a diagonal 
sub-block closed under commutation, isomorphic to a 
so(4) algebra. With O3 commuting with all the rest 
of this sub-block, the other six divide into two mu- 
tually commuting su{2) triplets upon forming appro- 
priate linear combinations which pass from a standard 
4x4 antisymmetric real-matrix representation of 5o(4) 
to 57t(2) ® 57i(2) matrix (4 x 4) form. Another such 
example, previously considered by one of us 0, is pro- 
vided by (O2, O3, Oi3, Oi4, Oi5, Oie), with O4 commut- 
ing with all of them. The linear combinations, defined 
in as "pseudo-spins", obey S7t(2) commutation rela- 
tions while not being isomorphic to Pauli spinors. While 
the Pauli matrices constitute a matrix algebra with the 
unit operator as the identity element, the pseudo-spin 
matrix algebra has as its identity element the operator 
^ ^/ ± (JzT^ . One implication is that this particular 
type of S7i(2) ® S7i(2) algebra is not related by a ST op- 
eration to the S7i(2) © 5it(2) algebras considered in the 
previous section. These 57t(2) © 57i(2) © 77(1) sets arise, 
for example, in the physical context of nuclear magnetic 
resonance |g with scalar coupling and heteronuclear two- 
spin coherences. For applications in methods such as 
unitary integration which require closure under commu- 
tation, such pseudo-spins work equally well as we will see 
below. There are many such sub-algebras which can be 
obtained by the discrete subgroup of ST operations as 
discussed in the previous section. For example, one can 
simply relabel {x,y,z). Indeed, each row (or column) of 
Table I has seven zero entries. The corresponding seven 
operators constitute such an algebra for a total of 15 sets 
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The recent demonstration of two qubit gates formed 
out of Josephson junctions |3| considers the Hamiltonian, 
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Expressing in terms of the Oi matrices in Table I, this 
Hamiltonian is a linear combination of (01,04,05,09). 
Together with (Ogj O12, 014), they form a closed sub- 
algebra, and the operator O13 commutes with all of them. 
For applications below in Section 2, we record the Hamil- 
tonian in terms of the Pauli spinors. 



H - ^iEoo+Ew)~^{Ej2ax+EjiT^)+^{Eoo-Ew)cr,T,, 

(4) 

and 



1 



{Eoo-E,o){s++s_-S+-S^ 



(5) 



where we have defined uj± ~ —Ej2 T Eji, — \{<Jx + 
Tx), Sz = \{(7x - Tj;), and 

S± = ]^{(Ty ± i(Jz){Ty ± iTz), s± = ^{ay ± icrz){Ty =F iTz)- 

(6) 

The sets (SztS±) and (5z,5±) mutually commute and 
each obeys the commutation relations of an 577(2) algebra 
but, as in they are not completely isomorphic to the 
4x4 matrix ai or algebras, as discussed above. In 
particular, as is easily seen, the squares of Sz and Sz are 
not the unit operator but involve well. While 

sharing the same commutation relations as the similar 
pseudo-spins defined in Q, the sets here differ as 4 x 4 
matrices, being in the "^-representation" instead of the 
counterpart "^-representation" in 



UNITARY INTEGRATION OF THE EVOLUTION 
OPERATOR 

Time-dependent quantum problems can be conve- 
niently handled through the unitary integration proce- 
dure of [ii . Solutions of the Schrodingcr equation for the 
quantum state are obtained through the evolution oper- 
ator U{t) defined by 



iU{t) = Hit)U{t), U{0) 



(7) 



where an over-dot denotes differentiation with respect 
to time. For more general treatment that includes 
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mixed states, master equations such as the Liuoville-von 
Neumann-Lindblad equation [l^ . 

= [H, p] + \^ {[LkP, Ll] + [Lfe, pL{]) 

k 

= [H. p] - \^ [lIUp + pilife - 2UpLl) , (8) 

k 

when first cast in terms of the individual density matrix 
elements, take a form similar to Eq. Q, namely 

Z77(i) - Ci^it), (9) 

where 77 are appropriate linear combinations of matrix 
elements of the density matrix pij ^J, such as, for exam- 
ple, the ones obtained by a decomposition of pij in the 
set of generators of the corresponding SU {n) group |0 . 
We have 




FIG. 2: As in Fig.l but for the real part of the off-diagonal 
elements of the density matrix. 



V = [pii - P22, (pii + P22 - 2p33)/V3, 

{pii + P22 + P33 - 3p44)/\/6, P12 -I- P21, 

i{pi2 - P2l),Pl3 + P31, i(Pl3 - P3l)- /9l4 + P41, 
i{piA - P4l),P23 + P32, i(P23 " P32), P24 + P42, 
«(P24 - P42),P34 + P43,«(P34 - P43)]- (10) 
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FIG. 1: Time evolution of the diagonal elements of the density 
matrix of a n = 4 two-qubit Josephson junction system. The 
horizontal axis is {ujt/2-K), parameters as given in the text 
with ui = 1 GHz, S — 0, and the initial configuration p(0) — 
SiiSji. 

In either case of Eq. {Tj) or Eq. lO, unitary integra- 
tion writes the solution as a product of exponentials, 
]^^- exp[— with time-dependent, classical func- 
tions Hjit), and Aj the elements of a complete algebra 
0. Upon substituting this form into Eq. ((TJ or Eq. iQ, 
a consistent set of defining equations for the fij follow 
so long as each element of the algebra is included in the 
above product. As shown in 0,13, so long as one has 
an su(2) algebra, solving this set defining pj reduces to 
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FIG. 3: As in Fig.l but for the imaginary part of the off- 
diagonal elements of the density matrix. 



just one non-trivial step, the solution of a single Riccati 
equation. This remains true if more than one su{2) is 
involved so long as they are mutually commuting, each 
set involving one Riccati equation. Note also that since 
only commutators of Aj are involved in the derivation, 
the results apply also to pseudo-spins as in the previous 
section, or for other algebras such as so(3) or su(l, 1) 
differing only in the sign or value of the structure co- 
efficients which do not change the form of the defining 
equations for the p. 

It is this feature that we now exploit in applying to 
four-level or two qubit problems that fall into the two 
categories of the su{2) © su{2) and su{2) © su{2) © u(l) 
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FIG. 4: Time evolution of the diagonal elements of the density 
matrix of a n = 4 two-qubit Josephson junction system. The 
horizontal axis is {u}t/2n), parameters as given in the text 
with uj = 1 GHz, 5 = 7r/4, and the initial configuration p(0) — 
SiiSji. Contrast with Fig. 1 to see differences due to a mutual 
phase between the driving fields. 
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FIG. 5: As in Fig. 4 but for the real part of the off-diagonal 
elements of the density matrix. 



sub- algebras discussed in Section 2. An example of the 
former is a recent discussion of geometric phases with 
two qubits. More general time dependences than con- 
sidered in that reference are amenable to our procedure 
presented below. For the coupled pair of Josephson junc- 
tions of |3| with the Hamiltonian in Eq. |(SJ , we have from 
our previous derivation in 5] the evolution operator 



U{t) — exp(— irio(i)) exp[— iij/_|_(i)S'+] cxp[— 
X exp[-iw3(t)S';,] exp[-ii^+(t)s+] 



2 

X exp[--i^_(i)s_] exp[--i//3(i)s^]. 



(11) 



The sets /i and ly are for the two SU {2)'s involved in 
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FIG. 6: As in Fig. 4 but for the imaginary part of the off- 
diagonal elements of the density matrix. 



problems such as those in An additional term in 
axTx in Eq. (01 could also be accommodated through 
an additional exponential factor in Eq. and would 
correspond to an U(l). These classical functions in the 
exponents arc obtained through the set of equations jj| , 



^'-f — /i+^/c- -|- iuj^ii^ 



(-E^oo + -£'io)/2, 

/is — 2ik^fi^ = Lu^, (12) 



with a similar set for the functions v which have K± 
instead of k± . These parameters for the H in Eq. jSJ take 
the values K± — —k± = (i?oo — -E'io)/2. For applications 
with k constant as in 3], and uj± functions of time, the 
first, and only non-trivial, equation of the set in Eq. \12l 
can also be recast as 



e{t) = k- \jj^{t)sm2e{t), 
with = tan(^^ + kt), or 

^{t)+lLO^{t)^{t)+k'-f{t)^Q, 



(13) 



(14) 



with = —{\/k){d/ dt){hi^). In any of these alterna- 
tive forms, (and, similarly, v^) can be solved for a 
given cj-t, following which the rest of the set in Eq. \i2\ 
yields to simple integration. 

In Figs. 1-3 , we present results for the two-qubit 
Josephson junction system of The parameters of 
that study were Eqq - Ew = 7.85 GHz, Eji = 13.4 
and Ej2 = 9.1 GHz. All these parameters were held 
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constant in that study but, with an eye to future ap- 
plications where these energies may vary with time, we 
consider such appHcations of our general time-dependent 
formalism. Specifically, were the voltages to be har- 
monically modulated, with Eji = 13.4cos(ijjt -I- 6) and 
Ej2 = 9.1 cos(Li;t), solution of Eq. ((T^ leads to diagonal 
and off-diagonal components of the density matrix shown 
in Figs. 1-3. Starting with an initial pure state (1, 0, 0, 0), 
the density matrix, p{t) = Up{0)W is evaluated and ex- 
hibits complicated oscillatory population of the four lev- 
els (diagonal density matrix elements in Fig. 1 which, of 
course, always sum to unity for all t) as also coherences 
(off-diagonal elements in Figs. 2 and 3). Figs. 4-6 show 
the effect of a phase difference between the two driving 
fields, leading to interesting differences when compared 
to Figs. 1-3 where the driving fields are in phase. 
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FIG. 8: As in Fig. 2 but with a damping constant F = 0.5 
GHz. Note the differences from Fig. 2 as the oscillations 
damp out asymptotically. 
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FIG. 7: As in Fig. 1 but with a damping constant F = 0.5 
GHz. Note the differences from Fig. 1 as the oscillations 
damp out to give asymptotic values of 1/4. 
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Extension to decoherence and the master equation in 
Eq. (jSl requires handling the 15 x 15 problem in Eq. Q 
rather than 4x4 and is necessarily more involved numer- 
ically. However, as in earlier studies 0,01 a simplified 
model in which the sum in Eq. JHJ) runs over all = Ok 
of Table 1 reduces to the same 4x4 calculation of unitary 
evolution above, with the only difference being to mul- 
tiply r\{t) in Eq. by exp(-Ft). Figs. 7-9 show the 
results of such a calculation for the same choice of pa- 
rameters as in Figs. 1-3 but with the additional F = 0.5 
GHz. The small t dependence of the two sets of figures 
coincide but the oscillations are damped by the decoher- 
ence such that asymptotically, the density matrix evolves 
to that of a mixed state, namely diagonal components of 
1/4 while all off-diagonal elements vanish. This evolution 
from a pure state to a completely mixed state is accom- 
panied by a monotonic rise of the entropy from to In 4, 



FIG. 9: As in Fig. 3 but with a damping constant F = 0.5 
GHz. Note the differences from Fig. 3 as the oscillations 
damp out asymptotically. 

paralleling the similar results in 0, 0| . 
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